from pyproj import Transformer

# 定义CGCS2000和WGS84的坐标系
# CGCS2000的EPSG代码是4490，WGS84的EPSG代码是4326
crs_cgcs2000 = "EPSG:4490"  # CGCS2000坐标系
crs_wgs84 = "EPSG:4326"     # WGS84坐标系

# 创建转换器
transformer = Transformer.from_crs(crs_cgcs2000, crs_wgs84)

# 示例CGCS2000经纬度点 (经度, 纬度)
cgcs2000_points = [
    (120.0839496198002, 68.19411326168881)
    # 添加更多点...
]

# 转换坐标点
wgs84_points = []
for lon, lat in cgcs2000_points:
    # 注意：transform方法的参数顺序是 (纬度, 经度)
    lat_wgs84, lon_wgs84 = transformer.transform(lat, lon)
    wgs84_points.append((lon_wgs84, lat_wgs84))

# 输出转换后的WGS84经纬度点
for i, (lon, lat) in enumerate(wgs84_points):
    print(f"Point {i+1}: 经度Longitude={lon}, 纬度Latitude={lat}")